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Abstract 

An asymptotic behaviour of solution of Kadomtsev-Petviashvili- 
2 equation is obtained as t — > oo uniformly with respect to spatial 
variables. 

1 Introduction 

Properties of solutions and process of integration for the Kadomtsev- 
Petviashvili equation (KP) M: 

d x (dtu + 6ud x u + d x u) = —3a 2 d 2 u (1) 

depend on sign of a 2 . This equation is called KP-2 if a 2 = 1. 

In this work the asymptotic behaviour of decreasing solution of 
equation KP-2 is obtained as t — > oo. The main term of the asymp- 
totics has an order by 0(t~ l ) and fast oscillates. An envelope of these 
oscillations depends on £ = x/t and r] = y/t. Results of the work 
are formulated in terms of scattering data for auxiliary linear prob- 
lem, which is associated with equation KP-2 in the inverse scattering 
transform (1ST) g |. 



*This work was partially supported by RFBR (97-01-459), grant for scientific schools 
(001596038) and INTAS (N0IO68). 
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The equation KP plays important role in modern mathematical 
physics. Many applications of these equation are known in plasma 
physics, water waves and other fields of waves process 0]. There- 
fore questions about solvability of this equation in different functional 
classes were studied in detail. In particular, the existence of global so- 
lution of Cauchy problem in class of distribution functions was proved 
in H|. A norm of solution in Sobolev space was estimated and an 
asymptotics of a solution were obtained as t — > oo for KP-like equa- 
tion but with more high nonlinearity, which called generalized equa- 
tion KP, in §]. 

The 1ST formalism allows to reduce a constructing of the solution 
of nonlinear integrable equation into solving of linear problems. One of 
most important achievements of 1ST is the constructing of asymptotic 
behaviours of solutions as t — > oo. These results are well-known for 
1+1-dimensional equations (one spatial and one temporal variable) 
@-|13]. Asymptotic behaviours of solutions for (2+l)-dimensional 



equations were studied not so in detail. Rigorous results about tem- 
poral asymptotics of solutions, which is nonuniform with respect to 
spatial variables, for a special class of nondecreasing solutions were ob- 



tained in [14 1, [15|. In work [16] an formal asymptotics of decreasing 
solution of equation KP-1, which is nonuniform with respect to spa- 
tial variables, was constructed. Results about asymptotic behaviour 
of decreasing solutions of equation KP-2 as t — > oo and about a re- 
mainder of this asymptotics, which are uniform with respect to spatial 
variables, are obtained in present work. 

The formalism 1ST for solving of Cauchy problem for equation 
KP-2 was presented in [fT^ [. This formalism will be used in this work. 
Below we remind basic steps of solving of Cauchy problem for equation 
KP-2 by 1ST. 

Let us denote an initial condition for equation KP-2 as: 

u\ t =o = u (x,y). (2) 

The process of solving of the Cauchy problem ([l|), @ consists of 
several steps. 

First step is solving direct scattering problem. On this step a 
boundary problem is solved for function <p: 

- dyip + dlip + 2ikd x ip + u<p = 0, y||jfc|-»ao = !> ( 3 ) 
and scattering data are constructed by a formula: 

F(k)= mt^mifj ■ 2d xdy UQ (x,yMx,y,k,0) x 
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x exp(-i(k + k)x - (k 2 - k 2 )y). (4) 



It is useful to note, if u is real, then the function F(k) has property 
F(—k) = —F(k). This follows from formulas (^) and (||). 

On the next step an evolution of scattering data is determined. 
This evolution is very simple: 

F{k; t) = F{k) exp(4it(fc 3 + k 3 )). 

Third step is solving of inverse scattering problem. This is reduced 
to so-called D-problem: 



djsp = ip F(—k) exp(itS), 
dk4> = —<pF(k)ex.p(—itS); 





(5) 



where S = 4(/c 3 + k 3 ) + (k + k)£ - i(k 2 - k 2 )r], £ = x/t, rj = y/t, the 
function F{k) is nonanalytic with respect to complex variable k G C. 
Solving of this problem allows to obtain the functions ip and ip at any 
time, if we know the evolution of scattering data. 

At last we can obtain the solution of the Cauchy problem by using 



the formula [17]: 



u(x,y,t)=d x J j^dk A dk F(k)ip(k, x,y,t) exp(itS). (6) 

However, this successive method of solving for Cauchy problems 
for general initial conditions allows to obtain very implicit answer 
with respect to initial data. Therefore the problem (|5|) and integral 
@ are used usually for solving of integrable nonlinear equations and 
the scattering data is taking from some suitable functional class (see, 
for example, [|l4|]-[]l6|]). On one-hand side this gives for studying a 
functional class of solutions, but on the other hand-side this leads to 
implicit restrictions on class of initial conditions. 

The auxiliary linear problem (H) and some related problems were 



studied in |T^]-[23]. In the work [^] it was proved that if the function 
uo(x,y) decreases exponentially, then the direct scattering problem 
df) is solvable for V7c G C. 

To construct the asymptotics of u as t — > oo we use the asymptotics 
of ijj and evaluate an asymptotics of the integral @ by stationary 
phase method As it turned out, the main term of asymptotics 

of the function u may be defined by using only the main term of 
an asymptotic expansion of ip as t — > oo. Taking into account the 
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fast oscillating coefficients of the system from (|5|) one can guess, that 
ip = 1 + o(l) as t — > 00. But we cannot say something definitely about 
the asymptotic behaviour of the function u, until we do not know 
analytic properties and an order of the remainder of asymptotics of ip 
more precisely. Thus we come to studying an asymptotic behaviour 
of a solution of the l)-problem (||). 

The asymptotic behaviour of the solution of the l)-problem as 
t — ► 00 with continuous and fast oscillated coefficients was obtained in 
1 25]. Here we study the Z)-problem with discontinuous coefficients on 
imaginary axis of complex parameter k. Constructing of the asymp- 
totics of solution for such problem is more complicated. Not only 
stationary points of phase function of an oscillated exponent define a 
structure of the asymptotic solution for (|5|) as in |25| (see also p6[), 
but the location of these stationary points with respect to line of dis- 
continuity of coefficients of the equations (||) as well. It contributes 
additional difficulties into evaluating and leads to changes in results. 
The uniform asymptotics of the solution of the problem (||) is con- 
structed by matching method [27]. 



2 Main result 



Theorem 1 Let (1 + \k\)F € L x n C, d a F e L x n C as Re(k) / 0, 
\a\ < 2 and: 

F(k + iX) 



sup 



dndX 



K + i\ — z 



< 2ir, 



(7) 



then the solution of the Cauchy problem of equation KP-2 for corre- 
sponding initial condition exists as Vt > 0. The asymptotic behaviour 
of the solution as t — > 00 differs in different domains of variables 
(x,y,t): 

as -(12^ + r? 2 )* 1 / 3 > 1; 



u(x,y,t) 



-At- 1 - 



7T 



x exp 

as (12^ + r? 2 )t 1/3 > 1: 



=f( 1 - 

Wit 



u 2 - x . 

^-12-) +c,-. + o(l). 



u = o(t -1 ); 
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as |12f + I2r] 2 \ <C 1; 
u(x,y,t) = 8it _1 0r/(ir//12)^ dpi cos ^8pl - -pi j + 



^ dpi^fpi sin ^8p? - zrH + o(t 



.Here £ = x/t, n = y/t, 



,( y 2 , x V 

\i2t 4 /3 + t i/3;' 



f( k ) = ^~ I l 9 dxdyu (x,y)ip(x,y,k,0)x 
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x exp(-z(/c + - (A; 2 - k 2 )y). 



The domains of validity for the asymptotics of the solution of equa- 
tion KP-2 are intersected and therefore they give combined asymp- 
totics of the solution uniformly on plane of x, y. 



3 An analytic behaviours of the scat- 
tering data 

In this section we demonstrate analytic behaviour of the scattering 
data which corresponds to sufficiently smooth and decreasing initial 
condition (|2|). 

First of all we show, that the solution of @ exists. 

Theorem 2 Let F(k) be such, that the condition (fy) is fulfilled, then 
the solution of (^) exists in a space of continuous vector- functions 
bounded when k G C. 

Proof. Let us consider a system of integral equations which is 
equivalent to (||): 

= ({)+G[F] 

where 

G[F]V = II dm A dm x 
J JmeC 

E«l „ m ^« , nm , ? ,„, t) . 





k—m 



exp(-itS) 



Using well-known results about integral operators (see, for example 
||29(|) one can show that the operator G[F] transforms the space of 
continuous vector-functions into itself. 

The operator G[F] is contracting operator. It is follows from in- 
equality (0). Hence we obtain the theorem statement. 

To construct the asymptotics we must study smoothness of scat- 
tering data F{k) in neighborhood of line of discontinuity Re(fe) = 0. 



Lemma 1 Let uq(x,v) be a finite function, then in neighborhood of 
some point k' £ C the scattering data have to be represented in the 
form: 



l\3\ 



F(k) = sgn[Re(-k)] £ f aia2 (k')(k - k') ai (k - fc')* 2 +0(\k-k' 
M=o 

where f ai ,a 2 (k') are continuous functions, when Re(k') ^ and 

W 2 (fc') = & 2 (k') + sgn[Re(-k)]fi% 2 (k'), 
when Re(k') = 0. 

Proof of this lemma consists of successive evaluating of partial deriva- 
tives of scattering data with respect to k and k. For example let us 
evaluate /{?(*/) and fH\k') as Re(k') = 0. 
Denote 

ThenfS ) (k') = f(k'),fS\k')^0. 

Evaluate f w (k') = f$\k!) + sgn[Re(- k)]f§\k>). 

dkf(k) = — f [ 2 dxdy u (x, y) exp(-i(k + k)x - (k 2 - k 2 )y) x 
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\(ix - 2ky)<p(x, y, k, 0) + d k <p(x, y, k, 0)] . 



To evaluate the derivative d k (p we use the integral equation for the 
function <p: 

8 W = ^V.P.[ f ^gF(-m)x 



2z7r J Jc (k — m) 2 
exp(i(m + m) — (m 2 — m 2 )y)'0(x, y, m). 
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The functions F(—m) and ip(x, y, m) are smooth on the left and right- 
hand sides of complex plane of m and one can show that the integral 
exists. Let us represent it into more convenient form. For this we 
rewrite the integral into sum of two integrals over left-hand side and 
right-hand side of complex plane. Let us integrate by parts these 
integrals. As a result we obtain 

**> = ^E/ [/<-*> * 

2m ^ Jdo,± k — m 
exp(i(m + m)x — (m 2 — fh 2 )y)ip(x, y, m)] — 
I (ix-2my)(±)f(-fh)ip(x,y,m) + 



2m ~^ J Jq± k — m 

(±)d m f(-m)ip(x, y, m) - tp(x, y, m)f(-m)f(m) x 
x exp(i(m + m)x — (m 2 — fh 2 )y). 

Here Q^ 1 is right-hand side (+) and left-hand side (-) of complex plane. 
We mean the integral over <9fi + as sum of integral over right-hand side 
and left-hand side of the circle with center at origin of coordinates and 
with large radius (which tends into infinity) and of improper integral 
over imaginary line and of an integral over the circle with small radius 
e with center at k = m. The direction of the integration over dVL + is 
determined by standard way. 

Consider the sum of the integrals over dQ.^. Each of the integrals 
over large half-circles is equal to zero as R — > oo because the integrand 
decreases. The sum of the integrals over the imaginary axis gives 
doubled integral over the imaginary axis in positive direction. The 
integral over small circle equals zero as e — > 0. 

Let us evaluate the derivative d m f(—m) of the integrand of the 
double integrals. For this we note that ip(x,y,k,0) = ijj(x,y,—k,0), 
hence /(— m) has to be represented by the function tp(x,y,m,0) and 
then: 

dmf{-m) = 7T- [ [ dxdyu (x,y)[- f(m)(p(x,y,m,0) + 



2vr 

+ (— ix — 2my)ip(x, y, m, 0)] exp(— i(m + m)x — (m 2 — m 2 )y). 
So we obtain a formula fio(k'): 

/io(fc') = fS(k') + sgn[Re(-fc)]/ 1 ( 2) (fc'), 
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where 



fil\k')=-f(-k')i;(x,y,k',0)^- , hr: 



dxdyuo(x, y) x 
{ix - 2ky)(p(x,y,k,0) + 

1 rice Afh, 



2ivr S / L 



dm A dm 



n± k — m 



{ix — 2my)(±)f(—m)ip(x,y,m) + 



(±)d m f(-m)ip(x,y,m) - tp(x,y,m)f(-m)f(m 

exp(i(m + m)x — (m 2 — fh 2 )y)^ x 
exp(-i(fc / + k')x - (k' 2 - k' 2 )y) 



fio = ~ I L2 dxd yuo(x, y)f(-k')ip(x, y, k', 0). 



An expression of foi{k') = foi(k') + sgn[Re(— AOl/oi^') has the 
form: 

foi\k') = 2^F / !^dxdyu {x,y){-ix + 2k'y)ip{x,y,k' ,0) x 

exp(i(fc' + fc'Jx - (fc /2 - fc /2 )y), 
fo?fi) = ^ I ls, 2 dxdyu {x,y)ip{x,y,k',0) x 

x exp(2i(fe' + k')x - 2(k' 2 - k' 2 )y). 

The expressions of other coefficients of the expansion have to be 
evaluated by the same way. The lemma is proved. 



f(2), 



4 Asymptotic solution of D-problem 

In this section we construct an asymptotic solution as t — > 00 of the 
Z)-problcm: 



dfr/j, = v F(k) exp(itS), ( n \ , 
d fe i/ = -//F(-ife)expHtS); I 1/ J H fc H°° 



1 





(8) 
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The solution of the problem (||) has to be obtained using the solu- 
tion of (JS|) and the formula: 

<p\ = ( fi(k,£,ri,t) \ ( -v(-k,£,r),t) \ , . 

The asymptotic solution of the problem (|8|) is combined. Here 
we construct the asymptotics as t — > oo uniformly with respect to 
all parameters. The stationary points of the functions S(k, k, £, rj) 
with respect to parameters k, k plays the important role in these 
constructions. 

The asymptotic expansion of the solutions is constructed using 
the inverse powers of large parameter t as a asymptotic sequence of 
asymptotic expansion outside of small neighborhoods of the stationary 
points of the function. Near the stationary points the asymptotic 
sequence has the form t~ n / 2 , where n = 0,1,2,.... Following the 
terminology of matching method p7f , we call the asymptotics outside 
of small neighborhoods of stationary points by outer expansion and we 
call by interior expansion the asymptotics near the stationary points . 
The right-hand side in the system of equations (||) is discontinuous on 
the line Re(fc) = 0, therefore the asymptotic solution is differentiable 
out of the line Re(fc) = 0. 

The domains of validity for interior and exterior asymptotic ex- 
pansions are intersected. This fact is used by matching method in 
order to construct unique combined asymptotic expansions. 

The phase function S depends on two parameters (£,77) £ K 2 . On 
the curve 12£ + rf = the confluence of two stationary points of the 
function S occurs. In this case we have one confluent stationary point. 
The structure of the asymptotic expansion of solution of (|8|) is changed 
here. As |12£ + r/ 2 | <C 1 the expansion is constructed on the powers of 
i~ n / 3 , n = 0, 1, 2, . . . as a asymptotic sequence. 

The formulas for the uniform asymptotic expansions of the solution 
of are large and it seems convenient to formulate the results about 
general case in section 4.1 and about confluence case in section 4.2 for 
convenience . 

4.1 Asymptotics in a general case 

In this section the combined asymptotic solution of (^|) is constructed 
when the phase function S has nondegenerate stationary points k = k\ 
and k = hi- Here we suppose that VA; G C, t l ^\df,S\k=k 12 \ ^ 1 ^ 
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t — > oo. This leads to restriction on values of the parameters £ and 
77, namely, i 1//3 |12£ + r/ 2 | S> 1. The asymptotic expansion which is 
uniform with respect to k £ C is formulated in the end of this section. 

To construct the combined asymptotic solution which is valuable as 
k £ C we obtain exterior and interior asymptotic expansions. These 
expansions are valid outside of small neighborhoods of kj, j = 1,2 
and in the small neighborhoods of kj respectively. Denote: 9 = 
\J — 12£ — r/ 2 . Then one can obtain an expressions for the stationary 
points of S by using the parameter 6: k x = ^(irj+O), &2 = 



4.1.1 The stationary points outside of the break line 

Consider the case when Re(#) 7^ 0, i.e. when the stationary points 
k\ t 2 are outside of the discontinuity line Re(fc) = 0. Let us formulate 
a result of this section about the combined asymptotic solution: 

Lemma 2 Let the system of the equations ^ have not the homoge- 
neous solutions, F(k) G C 2 D L\ as Re(k) 7^ and the parameters £ 
and r\ satisfy the inequality t 1 / 3 ^ 2 3> 1, then: 

when \/t\Q\\k — k\2\ 3> 1 the formal asymptotic solution of the 
system (^) with respect to mod(0(t~ 2 \dkS\~^)) has the form: 

fi = l + r 1 }i(k,$,r ) ), 

v = (r 1 v x (k, £, rj) + t~ 2 l x (k, £, 7])) exp(-iiS) + t- 1 u (k, f , 77); 
1121 



the functions V, v\, v x , u are defined by (fig), ^16j), 

when I A; — kj\ <C 1 the formal asymptotic solution of the 
system (^) with respect to moci(O(t _1 |0| _1 )) has the form: 

fi= 1 + M (lj,S,rj), 

v = (r 1/2 N (lj,t,v) N (lj,^,v))exp(-itS), 
where L ,j = 1,2, are defined by formula: 



lj = Vt(k-k j )J^ + 4(k-kj) 



112 
the functions M (lj,^,rj), N [lj,£,r]), N {lj,£,v) are defined by (\3_ 
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Proof. Let us construct the external asymptotic solution in the 
form: 

H ex = l + r 1 h (k, e, v) + t~ 2 Si (k, £, 77) exp(iiS) + . . . , (10) 



^ = (t" 1 ^ (k, e, r?) + t" 2 I (k, e, ??) + ... ) expHiS) + 

+t- l ^ (fc,e,r?) + .... (11) 

Let us substitute the formulas (|l0| ) and (11) into equate coeffi- 
cients with equal power of t. As a result we obtain: 

-id k S v x (k,C,ri) = -sgn(Re(k))f(-k), 
8- k H (k, £, v) + id~ k S li (k, £, r?) exp(z'iS) = 

-sgn(Re(fc))/(A:) £1 (£,£,77) + 
sgn(Re(/c))/(A;) ^ (k,£,rj)exp(itS), 
-id k S exp(-itS) v x (k, £, rj) + d k v (k, £, r?) = 

( - sgn(Re(fc))/(-fc) i (k,£, V )d k vx (k,^,rf)) exp(-itS). 

If we equate to zero the coefficients of oscillated terms and nonoscil- 
lated terms of the expansions ( JToj ) and (|ll|) respectively , then we 
obtain formulas: 

1 _ sgn(Re(fc))/(-fc) 



id k S 



(12) 



d k l (k,^,ri)=0; (13) 
d A= -^- mk) - (14) 

mi (^^) = sgn(Re(fe))/ i fe ^° (is) 



2 



9 S 5 



*= (-sgn(Re(fc))/(-*) J 5 fc uj) . (16) 

The function V\ has a jump on the imaginary axis of k. This jump 
is . In order for the coefficients of the asymptotics of function v as 
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t 1 to be continuous, we add an analytic function of k, which has the 
same jump on the imaginary axis in inverse direction, to v\ exp(iS't): 

a ck) = - r r dnf{ :l ■ (17) 

in J-ioo (k - n)d n S 

i i 
The function v \ define an analytic function vq of k uncompletely. 

The rest terms will be defined when we will match the exterior and 

interior expansion. 

A Cauchy-Green formula gives solution of (|14|) which is decreasing 

as \k\ — > oo and bounded when Vfc£ C: 

i -1 f f d P Adpf(-p)f(p) 

2m J Jc k - p id p S ' 1 ' 

We can obtain the domain of values of k, where the external ex- 

1 2 

pansion is valid, using the condition u\ j(t~ 2 v)\ 3> 1. As a result 
of calculations we obtain: 



t\0\\k - kj\ > 1. 

Let us construct the interior asymptotic solution which is valid in 
the neighborhood of point kj as j = 1, 2. Denote new scaling variable 
by If 

g = t(k-kj) 2 ^- + At{k-kif. (19) 
J 2 

When \lj\ is not so large (t 1//2 |(9f Sj| 3 / 2 » an asymptotic for- 
mula is valid as t — > oo 



{k -^ = im''-m^'^---- (20) 

Rewrite the system (||) into terms of new variables lj and lj. Sub- 
stitute the asymptotic expansion 

r = l+i _1 Mj (lj,Z,v) + ---, (21) 

v™ = (r 1 / 2 Nj (l j ,C,v)+t~ 1 Nj (lj^,v) + ■ ■ .)exp{-itS), (22) 
into system (||). 
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As a result one can obtain equations for the coefficients of expan- 
sions (21) and ([22]): 



d h Nj -2ih Nj= sgn(Re(fc,)) 



2 

16 



8 h Nj -2il 3 Nj= -^^-sgn(Re(%))/(-%); (23) 

2 



fio(-kj))lj 



21; 



foi(-kj) 



i 1/2 
^ Mj= -sgn(Re(%))/(/c i ) jVj W^r- 



(24) 



(25) 



Construct the solutions of the equations (23)-fl2q). Since the ex- 
ternal expansion has not the terms of order t^^> 2 , then the boundary 

l 

condition for function Ni (lj,S,,rj) has the form: 



0. 



The solution of the boundary condition for Nj (lj,x,y) is evaluated 
by formula: 



sis, 



dn A dn 
lj-n 



2m 



exp(-i(n 2 +n 2 )); (26) 



The solution of nonhomogeneous equation for the function Nj 
nas the form: 



Nf{lj,£,v) =sgn(Re(%)) 



8f(-kj) _ /io(-fcj) 

B w le " p(i(?+?)) 



2/oi(-A;j)exp(i(Z: ■ + /j)) /" y A 
: ^S~\ 2m 



lj-n 



■ exp(— i(n 2 + re 2 )) 



(27) 
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The partial solution of the equation for the function Mj (lj,£,v) 
has to be written as four-multiple integral: 



where 
J 



1 

2i^ 



-2f(k,)f{-k 3 ) 
dn A dn exp(i(n 2 + n 2 )) 



J, 



(28) 



u 



n 



2ztt 
dm A dm 



n — m 



exp(—i(m 2 + m 2 )). 



It is possible to reduce this four-multiple integral into two-multiple 
integral (see Appendix): 



J = 1 



exp(ig + m 



2iir 



dn A dn 



exp(— i(n 2 + n 2 )) 



exp(^ + Z|)). 



There exists an domain of values of complex parameter k, in which 
the external and internal asymptotic expansions of solution for the 
problem @ are valid. In this domain these asymptotic expansions are 
equal up to the terms of order o(i _1 ). In the domain, where i 1 / 2 !^ — 
kj\ 3> 1 and \k — kj\ < t~ 1//4 , the external and internal expansion are 
valid. We compute asymptotics of external expansion as k — > kj and 
internal expansion as \L\ — > oo. 



1 2 

Let us present the asymptotics of the functions Nj and Nj' 



as 



oo: 



Nj \M- 



sgn(Re(%))J^r/(-fc 



2iL 



+ 



k u 3 

exp(i(/ 2 + e)) 



u 



+ o(\i 



I — 3> 



Nj (lj,^v)\\lj\^oc = sgn(Re(fc,-)) 



8if(-kj) _ ifw{-kj) 
. (dlSj) 2 



+ 



2ih 



%8i 



1^1 
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The matching condition for v means that in the domain t 1//3 |fc — 
kj\ 3> 1 and \k — kj\ = o(l): 

(V 1 / 2 Nj (l J ,tv)+t~ 1 Nj (lj,S,v)j exp(-itS) - 
r 1 ^ (k,£,r])exv(-itS)+ u (k,£,Ti)) = oit- 1 ). 

2 

Using the asymptotics of interior expansions of Nj, which are rewrit- 
ten in the terms of external variable k as k — > fci2, and external 

expansion of Vq, we can obtain: 

b ° <*.*•»> -*(*) + ^ffgf + 

. -2/(-fc 2 )cxp(itS 2 ) /nq\ 

+ |8gS a |(k=fa) ' {ZJ) 
where the function vd(k) is defined by (FT 



2 2 



1 



Nj=N/ + ^ / (A; j )exp(i(Z| + Z|))exp(^5 i ) + 

, 2sgn(Re(fc m ))/(-Ai m ) exp(itS m ) 2 - 

+ io 2 g jTT , =7 exp(i(Z 3 - + l^.)), (30) 

2 

where m ^ j, the function Nj s is defined by 
l 

Let us construct Mj (lj,£,,r]). The asymptotics of (pq) as \lj\ — > oo 
has the form: 



\ 



dp, 12* dp 



i 

Evaluate the asymptotics of I^q as A; — > fcj-: 



1 ,j(-p)f(,p)-mm-k) 



here n£ 1,2, n / j. 
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The matching condition for /x when t 1//3 |& — kj\ 3> 1 and \k — kj\ = 
o(f -1 / 4 ) has the form: 

(1 + i" 1 {k, £, ??))(! + t~ l Mj v)) = o{t- 1 ). 

l 

These condition allows to define Mj (lj,£,r)): 

Mj rj) =M/(lj,t V) + Cj(Z, rj). (31) 

Here the function Mj s (lj,{,,rj) is defined by Cj(^,rj) has the 

form: 

1 /" y ^ A ^/(-P)/(P) - f(kj)f(-kj) 



2m J Jc 1 I2i(p-kj) 2 (p-k n ) ' V 7 

where nGl,2, j. 

Thus we have constructed the interior expansion in neighborhood 
of nondegenerate stationary point of S as |12£ + rj 2 \ > t -1 / 4 . 

Lemma is proved. 

Constructed asymptotic solutions are nonuniform with respect to 
k. But one can obtain an uniform asymptotic solution by using their 
combination. This uniform solution has the form (see, for example, 

w 



+ 



Mi ^ 


*( 


M 2 


N X j 


iV 2 




( Mi 
V JVi 





M*)- (33) 

The result of action of operator A n ^ on the function M is defined 
as follows (f2?J). Take the formula for M(lj,£,r],t), and change the 
dependence on variable lj into the dependence on variable k using the 
formula (19). Rewrite the sum of all terms of asymptotic expansion 
up to t with powers equal to — m, where < m < n. For example for 
the functions M(Zi,£, rj, t) and N(li,£,r],t) this process leads to the 
formulas: 



4 (Ms 1 . ,-i V ( k - fc 2 + 4 ( fc - fci) 3 itMm) 
lM lj + ^(k-hr + ^k-hr 2^5x1 ' 
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H-h) 

2ie^e{k-k l Y + ^k-ki) i 

\eye{k - ktf + ±{k - k^ 
\ _ /oi(-fci) ye(k - hy + 4(k - h)* 

" \e\ y/6(k - ktf + 4(k - h)* 



+ 



If we substitute (|33|) into (g) and evaluate a remainder, then we 
obtain: 

Theorem 3 The formal asymptotic solution of the problem ^) with 
respect to mo<i(0((i|6>|) _1 )), which is uniformly valuable when k £ C 
and 6 2 t > 1, has the form (&5 



4.1.2 The stationary points on the imaginary axis 

If Re(0) = 0, then the stationary point of the phase function S belongs 
to the line, where the coefficients of the equation are discontinuous. 
In this case constructing of the formal asymptotic solution of problem 
(||) differs from the asymptotic solution which was constructed before. 
The main result about the combined asymptotic solution is: 

Lemma 3 Let the system of equations (W have no homogeneous so- 
lutions, F{k) G C 2 fl L\ as Re{k) ^ 0; the parameters £ and r\ are 
-t 2 l z (Ylt, + rj 2 ) > I, then: 

when y/t\6\\k — k\$\ ^ lthe formal asymptotic solution of with 
respect to mod{0{t~ 2 \dkS\~ ?J )) has the form: 

fi = l + t~ l a (k, £, 77), 

v = (t- 1 fa (k, e, rj) + t~ 2 /?! (k, e, 7])) exp(-itS) + r 1 f3 (k, £, 77); 
1 1 2 1 

the functions a, (3 1 , /3 1 , /3 are defined by jpoj), fiSty, mW, §$0> 
when \6\ \k — kj\ -C 1 the formal asymptotic solution of QQ) with respect 
to mo(i(0(i _1 |#| _1 )) has the form: 

Y = i + r 1 Y (ij,£,v), 

Z = {t- 1 ' 2 Z (l^M+t- 1 Z (lj^, V ))exp(-itS), 

112 

where the functions Y (lj,^,rj), Z (lj,£,,n), Z £,rj) are defined by 
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The proof. Let us construct the asymptotics. The external ex- 
pansion is constructed similarly as in 4.1.1. The main difference is that 
the stationary points of the function S are on the line of discontinuity 
of the coefficients of system (^). It leads to sufficient modifications in 
the formulas for asymptotics. Let us find the external expansion in 
the form: 

fi ex = 1 + r 1 a (k, f , 77) + r 2 a x (t, £, r?) exp{itS) + ..., (34) 

v ™ = t- x {J} (k, e, 77) exp{-itS)+ /3„ (k, C, ??)) + r 2 I exp(-ziS) + . . . . 

(35) 

After substituting of (|34j) and (^) into the system (^) we equate the 
coefficients with identical powers of t and with oscillated and nonoscil- 
lated terms correspondingly. As a result we obtain equations for the 
coefficients of asymptotic expansions (p4l ) and (|35|). The obvious for- 
mulas are 

OL= — f [ dnAdn /(~ n )/( n ) C3QA 

2in J Jc k-n id n S(n,^,r]) 7 
1 _ sgn(Re(fc))/(-fc) 

1 

The function f3 1 is discontinuous on the imaginary axis. We add an 
analytic function with the same jump on the imaginary axis in back 
direction then the coefficient of the expansion (^) at t~ l is continuous: 

= 1 dnf(n)dnS = 1 dnfjn) 

nidfrS J-ioc (k — n)d n S nid^S J-ioc k — n 

1 

Unlike v'q, this function have singularities on the discontinuous line. 
The formulas for the other coefficients of the expansions ( |34"|) and 



(35) have the forms: 



2 sgn(Re(fc))/(fc) P Q (k,£,r,) 

ai= ' (39) 

^ = ^s(- s ^( Re (k)f(-k)»-dkk)- (40) 
1 

Here the analytic function (3 is still undefined. We will finally define 
this function after matching of external and internal expansions. 
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The interior expansion in the neighborhood of the point kj depend 

on the scaling variable lj. The asymptotics of the functions /i and v 

1 2 

has the same asymptotic sequence, but the equations for Zj, Zj and 
1 

Yj have discontinuous right-hand sides. 



d h Zj -2ilj Zj= , /^/(-^)sgn(Re([exp(m/4))); (41) 



2 2 - 

dn Zj —lilj Zj= sgn(Re[Z exp(i7r/4)]) x 



*(d 2 k s^ ho[ j) + d 2 s/ 0l[ 3) \d 2 k Sj\ 



+ w/^ { ~ h)+ \m fSi ~ k) - (42) 

d- h y,= sgn(Re(/exp(-^/4)))/(/c J ) Zj ^j^. (43) 

Continuous partial solutions of the equations (f41~|), ([42]) and (43) 
are obtained by Cauchy-Green formula. The formulas for the partial 
solutions of the equations (^l]) and ( |42| ) are differing from the formulas 
obtained in section 4.1.1 by function sgn(— Re(nexp(i7r))) under the 
sign of double integral: 

Zj s = J[Hj](lj), Zj s = J[Hj](lj), (44) 

1 2 

where Hj and Hj are the right-hand side of the equations (|41|) and 
(42). An operator J has the form: 

exp(i(^ + Tj)) r r dnAdn ^^^ , n2 t ^ 



Jim) = J J c W "(«) «phw + 9))- 



The partial solution of ( |43| ) contains of the double integral the 

same as in (^) and the term, which is defined by integral over the 

l 

discontinuous line. The function Yf has the form: 

Y! = 2/( ,^ fcj V " 2J + - - 2J_ + ), (45) 



where: 



, . dm A dm exp(i(m 2 + m 2 )) 
' /() ■ I — m 1m 

dnAdn 2 _ 2 

exp(— i(n +n )). 

n± n - m 
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Here fi± = {±Re[Z exp(-w/4) > 0}. Note that J-+(l) = J+-(-l)- 

In the Appendix we show, that we can reduce the integral J |_ into 

double integral. As a result we obtain: 



J- 



J_ 



3. 

--Z7T + 
4 



il 



dm A dm 



n+ il — m 



exp(—i(m 2 + m 2 )), when I G 



Zexp(i(Z 2 + P)) 



e?m A dm 
n+ il — m 



exp(— i (m 2 + m 2 )) 



-—iir——i7rexp(i(l 2 + P)), when I G 17 ; 



Thus, the function Yj 5 is represented by sum of double integrals. 

Let us consider matching of the external and internal asymptotic 
expansions. For that we need the asymptotic behaviour of the partial 
solutions of (fTI|) and ( f42|) as \l\ — > oo. The asymptotic behaviour of 
the double integrals is evaluated in the Appendix. By using it we 
obtain: 



Z/\\i\- 



2 

Z 3%\~ 



sgn[Re[-^ exp(m/4)]]f(-kj) x j ( — + 



■1 



+ 



exp(-i(l? + i?)) 



21; 



+ OQL 



sgn[Re[— lj exp(i7r/4)]] x 



+ 



+ 



i i 

4 ~ 4 



exp(z(Zf + Z|)))2/ ( 



(i) 

01 



+ 



\%Sj\ 



+ exp( , (/ 2 + ^ + i/( g) + (| Zi |-i). 



The asymptotic behaviour of the function /3q' is: 



/5q |fc— >fcj 



+- 



/(tj)5gn[Re(-*)] + 



-oo n — ikj 



+ 
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+/ ( g ) (fc i )sgn[Im(iA ; )] + fio(kj) + /^(fc>gn[Im(ifc)] + 



— : / dn 



V2i{kj kff 
f(in) ~ f(kj) 



/(%)sgn[Re(-fc)] + 



n — ikj 



+ 



1 



(2), 



/o ( i 1} (%) + 



+ - 



7T/ 



00 dn /(w) - /(fcj)(m - ik^fjkj) - (in - ik^f^jkj) 
-oo (n — ifc,) 2 

+o(l). 



c(2) ( 



Using the matching conditions for the asymptotics of the external 
z/ ea: and internal v tn expansions we obtain: 



/Wo 



12*(fc - fci)(A; 1 - k 2 ) 12»(jfe - k 2 ){k 2 - h) 



(46) 



where (5q is defined by formula 



C} = - 
l l 



-oo n — 



Zj=Z S j -exp(-i(l] + tf))C], 



1 2 . . 

where functions Z/ and Z/ are defined by ([44J 



(47) 



1 ^ _J_ f (2) 



/lO ( %) 



1 

7T?' 



r(2)/ 



(n — ifcj) z 



The asymptotic behaviour of the function Y f is : 



1 

Yf 



m + 



\f(-kj)f(kj) 



12vrL 



The matching condition for the external expansion of \i ex and for the 
internal expansion of fj, tn gives: 



1 1 



./("%)/(%) 



(48) 
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1 

where Yf is denned by (f45|), The function Cj{^,rj) is defined by (|32|) . 
The lemma is proved. 

Constructed external and internal expansions are irregular with 
respect to parameters k, £, rj. The uniform expansion with respect to 
k when(12£ + r/ 2 )^ 1 / 3 S> 1 has to be constructed by the same way as 
in preceding section: 

(49) 

The operator A n ^ was defined above. Let us substitute ( |49| ) into (||) 
and evaluate a remainder. As a result we obtain 

Theorem 4 The formal asymptotic solution of the problem ^) with 

respect to mod{0({t\6\)~ 1 )) , which is uniform when k £ C and 
_ 6/ 2 f 2/3 ^ : has the j orm 



4.2 Asymptotics in neighborhood of confluent 
stationary point 

The system of equations in the problem (||) depends on two control 
parameters ^ and rj. On the parabola 12£ + rj 1 = the degeneracy 
of stationary points occurs: k\ = k<i = k$ = in. For this reason 
the asymptotics constructed in section 4.1 is invalid when parameter 

6 = \/n 2 + 12£ is close to zero. For example, the asymptotics of I 1 as 
k — > ko and 9 — > is discontinuous: 

■i -I / | k — ko | 9 k — ko 

A* Ifc-fcoJ |e^o = ^ 12i (( fc _ jfc )2_^2) " ^i2i((A;-A:o) 2 -0 2 ) 

f(-h)f(k ) + 

This shows that we need new scaling of the parameter 9. The scaled 
control parameter is: 

v = t 1 / 3 

The external expansion constructed in sec. 4.1 becomes discon- 
tinuous at 9 = and the internal expansions constructed in sec. 4.1 
become singular at the point 9 = and lose their asymptotic prop- 
erties. Therefore here we must change the internal variable for the 
internal asymptotic expansion: 

p = t^ 3 {k - k ). (50) 



V12 
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In this section we construct a formal asymptotic solution of the 
problem (H) with respect to mod{0(t~ 1 )) when \0\ <C 1 uniform with 
respect to k £ C. The result is formulated in the end of this section. 

To construct the uniform asymptotic solution we need the external 
and internal asymptotics outside and inside a small neighborhood of 
the point ko respectively. 

Lemma 4 Let the system of the equations (^)be have no the homoge- 
neous nontrivial solutions, F(k) £ C 2 D L\ and the parameters £ and 
r) satisfy the inequality |12£ + n 2 \ <C 1, then: 

when Ik — kolt 1 / 3 3> 1 the formal asymptotic solution of the system (^) 
with respect to mod(0(t~ 2 / 3 j\k — k§\) + 0{tr v )) has the form: 

fi = l + r l m (k,£,rj), 

v = t~ 2 / 3 n +t~ l ( rii ex.p(itS)+ n ) 

the function m is defined by ^/h]), the functions n and n-2 are defined 

by ^T^j and ( f7^ j respectively, the function hi is defined by $6^); 
when \k — ko\ <C 1 the asymptotic solution of the system ^)with respect 
to mod(0{t~ 2 / 3 \k - k \ + O^ 1 )) has the form: 



l + t- 2/3 M +t~ l M, 



v = (r 1 / 3 jtf +t~ 2,Z N +t- 1 AT) exp(-itS); 



the functions Jv[, j = 1,2 and J\f, j = 1,2,3 are defined by (75), and 
(|7^j respectively. 



Proof. Let us find the internal formal asymptotic expansion for 
the solution of the system (g) in the form: 

M in = l + t- 2 ' 3 M +t~ x M +..., (51) 

M in = (t~ 1/3 M +t" 2/3 M +t~ l Af + ...) exp(-itS). (52) 

Change the variable k into the variable p in the system (^). The 
phase function of the exponent is depend on the variable p as: 

tS = uj(p) = 4(p 3 +p 3 ) - v 2 (p + p), 
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where v = t^l^O / \/V2. Substitute (^Tj) and ( |52"| ) into the system (||), 
equate the coefficients with equal power of t. As a result we obtain 
a sequence of equation for the coefficients of the expansions (|5l]) and 



(52). 



d p AA -i(12p 2 - v 2 ) M= sgn[Re(-p)]/(-fc ), 
8 P M -i(Ylp 2 - v 2 ) M= -sgn[-Re(p)](f^(-~k )p + 
+A 1 o\-~k)p) ~ f$(-ko)P ~ f$(-ko)p. 
d p M= sgn[-Re(p)]/(*o) J\f, 
d p TV" -i(12p 2 - v 2 ) M= 



(53) 

(54) 
(55) 



sgn[ - Re (-p)]/«(-fc ) + f$(-h 



(2), 



sgn[-Re(-p)]/«(-fco) + / ( S2 ) (-^o) 



sgn[ - Re(p)]f[{\-k ) + fg(-ko) 



\p\ 



-sgn[ 

2 



Re(-p)]/(-*„) M 



(56) 



d p M= sgn[Re(p)]f(k ) M + 



sgn[-Re(p)]fH\k ) + f^ 



e(2), 



+ 



sgn[-Re(A ; )]/ ( S 1 1) (A ; o) + /^ 



f(2), 



l 



(57) 



These equations are obtained in a supposition that |p|t -1 / 3 <C 
1. The uniform bounded partial solutions of the equations for the 
coefficients of the expansion (|52] ) are obtained by using an integral 
operator: 



2iir 



dr A dr 



p — r 



exp(-iu(r))g(r), 



(58) 



where g(r) is the right-hand side of corresponding equation. 

Thus, the bounded partial solutions of the equations ( |53[ ) and (54) 

1 2 

(the functions J\f s and J\f s ) are represented by double integrals. 

Bounded partial solutions of equations for the coefficients of the 



expansion (|5l|) are constructed by using the Cauchy-Green formula. 
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The formula for the partial solution of the equation ([55]) has the 
form: 



M s (p,d,v) = f(k )f(-~k )(Ji(p,v 2 ) -2Jf+ -2J+-). (59) 

Here 

1 /" /" dn A dn exp(iu;(n)) f f dr A df 

Jl = / / / / pyp(— ?.f.;(r)). 

2in J JC P - n 2iir J J£ n-r 

T=c+ Iff dn A dn exp(iw(n)) /" f dr A df , . 

Jr = / / / / expf— aa;(r)). 

1 2 ? W p-n 2ztt J Jn± n-r ^ y " 

In the Appendix the integrals Ji and J-f + are reduced into double 

integrals. By the similar way we may represent the integral as 

l 

the sum of double integrals. Therefore we may represent J\A S as sum of 

double integrals also. Corresponded formula is very large and we don't 

write it here. But we will use this formula to evaluate an asymptotic 
l 

behaviour of M s as \p\ — > oo. 

The bounded partial solution of the equation (|57| ) when p£C has 

2 

to be built similarly . We give following statement about M- 

2 

Lemma 5 Continuous partial solution of the equation for M. exists. 
This solution is uniformly bounded when p£C. 

2 

Sketch of proof. The partial solution of the equation for M. has 
to be obtained as a result of using of the operator: 

„ i , , 2 1 f f dn A dn , . . . 

dp 1 [g] =M S ^— I g(n) 60 

y 2z7r J JC p — n 

to the right-hand side of the equation. The integrals, which are ob- 
tained, are continuous by virtue of continuity of the integral operator 
with respect to the parameter p. The boundedness of these integrals 
with respect to p when p € C may be obtained by using the asymp- 

2 

totic behaviour of the right-hand side terms of the equation for M as 
\p\ — ► oo. 

Let us to construct the external expansion in the form: 

1 2 2 

m = 1 + t' 1 m (k, v) + r 5/3 (mi exp(itS)+ m ) + (61) 
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n = (t" 1 hi (k, v) + t~ 5 / 3 rii (k, v) + t~ 2 n x . . .) exp(itS) 

+r 2 / 3 h (k, v) + r 1 h (k,v) + (62) 

Substitute ( pl| ) and ( |62"D into @. As a result we get the equations 
for coefficients of the expansions: 

A 1=S gn[Re(-fc)1 .ff"*? , 2 ; n 1= - S gn[Re(-fc)] ''^^ • 



12z(fc - A: ) 2 ' J1 lUi(k- & ) 4 ' 

(63) 

~ i f(k)f(~k) 

d * m = m(k-k r (64) 

5 fc n (fc, C, r?) = 0; <9 fc n (fc, £, r?) = 0; (65) 



_ 2 f(k)f(-k)v 2 

^ m ° = 144z(fc-fc )4 ; ^ 
% = sgn[Re(-fc)]/(fc) i ^- ( ^^ ) ) 2 ; (67) 



3 I - - 1 

(fc,f,f7) = 12j / fc _ fco p I sgn[Re(-fc)]/(-A:) m 
2/( - S) .^- t) ]j^zBg^(z^g(zBV (08) 



12i(fc - k f 12*(fc - fc ) 2 

12 3 

The formulas (|6^) and (|68| ) define the functions rti, n x and n x . Us- 

1 2 

ing the formulas (|65| ) we can see, that the functions uq hq are analytic 
of variable k. The obvious form of this dependence is defined by two 
conditions. First one is the continuity of the asymptotics and second 
one is matching condition for the external and internal asymptotic 
expansions. 

One can see, that the sufficient condition for the continuity of 
the coefficient of the asymptotics ( |62"D as t~ 1 with respect to k is an 

2 

addition into no of the term: 

n = — 2 / T=^=/H)- ( 69 ) 

m 12i(k - k ) J ~i°° k ~ X 
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It follows from the formulas ( p7| ) , that the coefficient of the asymp- 

totics ( |6l| ) as i~ 5//3 is defined after evaluating of no, i.e. after matching 
of the coefficients as i~ 2 / 3 of the external and internal expansions. 

Consider the problem for defining of m in detail. This function 
satisfies the boundary condition: 



m o \\k\- 



0. 



The external solution is valid outside of small neighborhoods of the 
points kj, j = 0,1,2. Therefore solutions of the equation (64), which 
are smooth, bounded and decreased as \k\ — > oo: 



dr Adf f(-f)f(r) - f(-k )f(k ) 



C k 



Yliir - k n ) 2 



+ 



+ 



(k - k )f(k )f(-k ) 
12i(k - k ) 2 



(70) 



are defined to within an analytic function with respect to variable k, 
which has poles when t l ^\k — k$\ -C 1. The full definition of m will 
be done by matching of the external and internal expansions. 

For matching process we need asymptotic behaviour of the partial 
solutions of equations (|53l)- (|56|) as \p\ — > oo. Constructing of these 
asymptotics is reduced to evaluating of an integrals with weak sin- 
gularity in the integrand. Evaluating of these integrals is done in the 
Appendix. Here we write the asymptotic behaviour of partial solutions 
of equations (|53|)-(|56|) as \p\ 



oo. 



AT S (P,^V) = ^exp(^(p))/(-fco)[0oo(^ 2 ) " </>o> 2 )] + 

+± I eMMp))f(-~kMi(v 2 ) ~ 4> m {v 2 )] + 

+^/(-Msgn[Re(- p -)](l - 6XP( ^ (P)) ) + 0(\v 2 \\ P n. 
Hi p z p z 



P 



1 exp {iui (p) ) \ sgn [Re ( —p) ] 
p P ) 12i 

\ exp(iuj(p))(j)io(v 2 ) + 
p \lip 



fil\-ko) + 



/£>( 



-ko) + 
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expMp)) _ ^ 2)) + (p__ exp( W (p)) 



P 



P 



sgn[Re(-p)] 



12* 



P 



\2ip i 



f(l\- ko ) + o(\ P r)+o(\vf\p\- 



p 



: exp(io;(p))0oi(f 2 ) + 



M s (p,Z,v) = -f(h)f(-ko) 
1 
P 



-p 



+ 



I2ip 2 

12^ i(v 2 )^p) +iv 2 (i )00 (v 2 )<j )00 (v 2 ) ) - 



+ 12i[^ 1 (v 2 )^l(v 2 ) + 0oi(^)V-O + l(^ 2 ) - oV'Ol 



+ 



0(\v\ 2 \p\- 2 )- 



M s (p,Z,v) = f(ko)f$(-ko) 
f(k )fH\-k ) 

/(fco)/ (2) (-/c ) 



p 

-t— + 

\2ip 




( P + 

V Uip 




m 


P 2 




24ip 2 


2Ai 



sgn[Re(p)] + 



f( ko )fU(-k )J^ + <S>(v 2 ) + 0{\v\ 2 \p\~ l ), 



N s {p,i,v) 

sgn[Re(-p)]( 



1 



sgn[Re(-p)](l - exp(iuj(p)))fH\-k ) 



24i 

p 2 exp(iuj(p)) (i) 



12ip 2 



24i 



)/#(-**>) + 



p 2 A 2) (-feo) 

24ip 2 



+ 

+ 
+ 



/o(-^o)/(A ; o)exp(^(p))M/(t, 2 ) + OMpl' 1 ). 

The functions <E>(t> 2 ) and ^(v 2 ) are smooth and uniformly bounded 
when v 2 6 R. Here we use notations: 



tRcW>o 



dr A drr m r n exp(io;(r)), 
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dr A drr m r n exp(— iu(r)), 



=FRe(r)>0 

<Pmn= J J^dr Adfr m f n exp(iuj(r)), 

ipmn = J J^dr A dfr m f n exp(—iu(r)). (71) 

Evaluate an asymptotic behaviours of the coefficients of the exter- 
nal expansion as k — > ko and df,S = o(l). 



1 

n i \k^k 
1 



/(-fco) 
12i(& - k 

(i) 



12i(/c - k 

,(2), T x ,(2), 7 y fc-fcp 

Jio l _fc oJ + /oi (- fc oJ fc _ fc 



2Sgn[Re(-fc)] + 
/iV(-fc) + /S ) (-/-o) 



k — ko 
k — ko 



sgn[Re(-fc)] 



(!)« 



k - k 1 ^,(1) (k - k f 



k — ko 



02 



(k - kof 



+ 



+ 



+ 



^/ 2 ( ?(-^o) + / 1 ( ?(-M^ + 



jA 2) (-^o) 



(fc - fcp) 2 

(k - ko) 2 



sgn[Re(-fc)] + o(l). 



2 / 

n 



24vr(fc - fe ) 2 



7Ti/(-fc )sgn[Re(-p)] + V.P. 



f 



d\f(i\) 
A — iko 



+ 



12tt(/c - fe ) 



TTi/o^C-fcd) + TTiffti-ko) + (7Ti/l (2)(-fc0) + 



+7ri/(S 2) (-feo))sgn[Re(-fco)] + 



+V.P. 



dA/(iA) 



+ 



l 

12? 



vri/io } (-fco) + 



^ (A-ifeo) 2 - 

+ 27Tz/ 1 ( 1 1) (-fc ) + TTif^(-ko) + (7Tz/ 2 ( ?(-fc ) + 2nif[f (-~k ) + 

+^ 2 2) (-fco)) S gn[Re(-fe)]+y.P. f° dA/(a) 

J — oo 



(A - ifco) 3 



+ 0(1). 



- s -iW^F /( -' o)/(fco) 
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i2,(.-. )L (/ - ( -^ )sgn[Re( ^ )] + 

+fil\-~ko))f(ko) + f(-~k )(fH\k )sgn[Re(k)] + f$(k )) 



f(koKf^\-k ) S gn[Re(-k)] + (-**,)) + 



(2), 



12i(fc - k ) 2 

+f(-~k )(f^(k )sgn[Re(k)] + f$(k )) 



+ o(l). 



Let us do matching of the external and internal asymptotic ex- 
pansions of the function v. The matching condition for v in domain 
t -1 / 3 <g | A; - k \ < 1 when |0| < 1 has the form: 

(t~ 2/3 n n Tii exp(z'iS')) ~ 

Let us equate the coefficients with equal powers of the large parameter 
t. As a result we obtain: 



l 

n = 



1 



k — ko 



f(-koM (v 2 )-<t>oo(v% 



(72) 



2 

n = 



i 



fc— fco 

1 fCl)/ ,[.,+ , „2i,r / ,.2 



,_^(-*o)[^(^(^)] + ,_, ii 

1 f(l) 



/io 

+ ^^ ) (-fco)0oi(« 2 ), 
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«o)0ioKJ + 



/o\ ; (-^o)[^ + iK)-^oiK)] + 



where 0, ?/>, V are defined by ([71]); 

AA= V[g 3 ]f 2 (-~k )f(k ) exp(iuj( P ))^(v 2 ), 



(73) 



(74) 



where gj are the right-hand sides of the equations for J\f, the operator 
V\g\ is defined by fl5q), the function ^(u 2 ) is smooth and uniformly 
bounded with respect to t> 2 when v 2 £ M. 

Matching condition for the function /2 in the domain t -1 / 3 <C \k — 
fco | <C 1 when |#| <C 1 has the form: 

(1 + m ) - (1 + r 2/3 M ) = oir 1 ). 
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As a result we obtain: 



l 1 



M=M S ; M=M S n, 



(75) 



where M s is denned by (59), the function M s is defined by ( |60| ) and 
the function <&(v 2 ) is smooth and uniformly bounded when v 2 £ K. 



1 1 s 
m=m 



1 



k — ko 
2 



12i0 O i(f 2 )</>io^ 2 ) + iv £ <j> 00 (v 2 )(f>oo{v 2 ) 



k hi 



12iU+ K)^oiK) + </>oiKWiK 



01 



(76) 



Here m s is defined by (|70|) and the functions ■0 =t > ^ > 4> an d ^ are 
defined by formulas (fn|). 

Thus we have matched internal and external expansions of jl and 
v. The lemma is proved. 

These expansions are ununiform with respect to k. Now we can 
construct uniform asymptotic expansion with respect to k G C. Fol- 
lowing the matching method the uniform expansion is: 



ft 

v 



m 
ft 



+ 



Mi 
Mi 



A 



i,fc 



Mi 



(77) 



Here the operator A n ^ processes on the function Mi in the formula 
for Mi by followed manner. One must change the variable p into the 
variable k using the formula (^) and write all terms of the asymptotic 
expansion with respect to t with the powers are equal to — m, where 
< m < n. For example, one can obtain for the function M(p, £, rj, t): 



A lik [M(p,Z,ri,t)] = l + t 
1 



/(W(-*o) 



-k — ki 







k — ko 
2 



12i(k - k ) 2 
12i4> i(v 2 )4>io{v 2 ) + iv 2 (j)oo(v 2 )(f)o (v 2 ) 



+ 



k _ k ^^oo^ 2 )^^ 2 ) +^oo(^)^ + o(^ 2 
+I2i^ l {v 2 )^ l {v 2 ) + ^oi(^ 2 X(^ 2 )) " \^ 

^)^-^iS) 



+ 



01 
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•2 



+/(W^ , (-fc)( S ^ji-Si)'HW*)l). 

The obviously formula for ^^[AA] is more large and doesn't shown 
here. 

Theorem 5 The formula fir\ ) gives the asymptotics solution of the 
problem ^) with respect to mod(0(t~ 1 )) as t —* oo. This asymptotic 
solution is uniform with respect to k G C and \9\ ^> 1. 

5 Justification of the asymptotics of 
the solution of Z)-problem 

In this section we prove that the remainder of the asymptotics has 
order by i -4 / 3 uniformly with respect to k G C and this remainder has 
to be differentiable with respect to x. We call by the remainder of the 
asymptotics the difference between solution of the problem(p|) and 
constructed asymptotic solutions (|33|) when 9 2 t~ 2 ^ 3 S> 1, (|49"|) when 
_Q2 t -2/3 > 1 and when ^1 < l The differentiability of the 

remainder will be important when we will construct an asymptotic 
behaviour of solution of the equation KP-2. 

Theorem 6 Let d a f(k,k) G L\ n C 1 when \a\ < 2 w/ien is out of 
the imaginary axis and 



sup 



dk A dk , „ , , . . 
|F(fe | 



|/c — z 

then the solution of the problem is: 



< 2tt, 



: = s 



+ 0(W 3 ), (78) 



lo/ien G C, £, 77 G R. T/ie remainder of the asymptotics has to be 
differentiable with respect to x. 
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The proof. Let us write the system of differential equation for 



the remainder. Denote the remainder in ( |82[ ) by V. Substitute (82) 
into (||). As a result we obtain: 

( d- k 0\ / F(-k)exp(itS)\ 

{ d k ) V ~ \ F(k)eM-itS) ) V + L (?9) 

^lifc^oo = 0. (80) 

We denote by the vector / the residual which originates in (§) when 
we substitute the column into this equation (/2, v) T : 

fi = —dyl + F(—k) ex.p(itS)i>, 

f2 = —d^v + F(k) exp(—itS)p,. 

Let us denote by X the space of bounded and continuous with 
respect to k vector-functions with the norm: 

\\W\\ = sup \W\\ + sup | Wb| - 

feeC,(C,r?)eR 2 keCM,v)eR 2 

Consider a system of integral equations instead of the problem 

V = G[F]V + fl", (81) 
where G[F] is the integral operator: 



G[F)V= f I dm 



A dm x 

eC 



^exp(^) ( 

|MexpHt5) ) S ('"-'/■/). 

H= / dm A dm | /2 fo^ |f) 

J Jm£ l L \ , 

\ k—m 

Using obvious forms of the functions fi t 2 one can proof, that 

\\H\\ = 0(t~^ 3 ) 

uniformly with respect to £, rj £ K when € C. 

The operator G[.F] is contracting in the space X, therefore the 
solution of the integral equation (^l]) exists in X and is evaluated by 
0(t~ 4//3 ) uniformly with respect to £,77 £ R. 
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Show that the remainder of the asymptotics is differentiable on 
x = i£. Differentiate with respect to x the system of the equation for 
the remainder. Denote the derivative of the vector V by x- Then we 
obtain: 

X = G[F] x + d x G[F]V + d x H. 

The terms d x G[F]V + d x H may be evaluated by order 0(t~ 4//3 ). 
The operator G[F] is contracting, therefore one can obtain: ||x|| = 
0(t -4 / 3 ). The theorem is proved. 



6 Solution of the equation KP-2 

An asymptotics of the solution of the problem (|5|) as t — > oo may be 
written as: 

(t) = {*) {k ^^ t)+ {7) ~ k ' e ' * t] + or4/3) - 

(82) 

The second term in this formula is the solution of the problem (||) 
with the boundary condition: 





\k\— >oo 

The proof of the theorem 1. Let us substitute the function 
ip{k, k, £, rj, t) into the formula for the solution of the equation KP- 
2(^). Differentiate the integrand with respect to x. The main terms 
of the integrand are the terms which appear after differentiating of the 
exponent. The derivatives with respect to x of others factors of inte- 
grand are small because they depend on x slowly. The integrals of such 
terms over the plane are evaluated by the order 0(£ -4//3 ) uniformly 
with respect to (x, y) E K. Rewrite the main term of the integral as the 
integral over real plane: (k, A) G M 2 , where k = Re(/c), A = Im(fc). 
As a result we obtain: 

u(x,y,t) = —A I l dK A dX \k\/(k -MA) x 
J JM 2 

x exp(it(8K 3 - 24kA 2 + 2< + 4kAt?)) + 0(t~^ 3 ). (83) 

Thus the main term of the asymptotics of the solution of the 
Cauchy problem for the equation KP-2 is given by integral with fast 
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oscillating exponent. Let us evaluate the asymptotic behaviour of this 
integral. 

Let rj 2 + 12£ > and t 1//3 |r/ 2 + 12£| » 1, then the stationary points 
of the exponent are: = ^±\/rj 2 + 12£, rei^ = 0. The stationary 
phase method (see for example [|4|]) gives: 

u(x,y,t) = o(^ 1 ). 

Let rj 2 + 12f < and i 1 / 3 ^ 2 + 12£| > 1, then A 1)2 = ife, Ki,a = 
±i\/— r/ 2 — 12£. The asymptotic behaviour of the integral is: 

-i 1 f./i./_„2 , *1 



exp ^ - llitJ-^ - 12 + c.c. + o(i -1 } 



To evaluate the main term of the asymptotics of the integral when 
|12£+77 2 | = o(l), we substitute the scaled variables p\ = i 1//3 Re(/c— ko), 
P2 = t 1//3 Im(A; — ko) and the parameter t> 2 = t 2 / 3 (r] 2 + 12£)/\/l2 into 
integral (|83|). As a result we obtain: 



u(£,rj,t) =4it 1 f(k ) J J^ 2 dpidp 2 x 

pi exp(i(8p? - 2u 2 pi - 24pyji)) + 

Let us integrate the internal integral with respect to the parameter 
P2, use the even property of the integrand with respect to p\. As a 
result we obtain: 

u(x,y,t) = m- x ^f {iril 12 )(^ d Pl ^pl cos (8p 1 (p 2 l -8v 2 )) + 

+ ^ °^iv^rs in (8Pi(P?-8^ 2 ))) + o^ 1 ). 
The theorem 1 is proved. 

A Asymptotic behaviour of double in- 
tegral with weak singularity of the in- 
tegrand 

Here we obtain the asymptotic behaviour of integrals which are ap- 
peared when the asymptotic solution of (||) was studied. 
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Evaluating of the asymptotic behaviour of one-dimensional inte- 
grals with weak singular integrand and fast oscillated exponent was 



done in [24] p. 26 and [Eg] p. 332. The asymptotic behaviour of many- 
dimensional integrals with fast oscillating exponent was studied in 
|24||, |pQ|. The asymptotic behaviour of the Cauchy integrals with fast 



oscillating exponent in one-dimensional case was studied in [ 24 1 and in 
many-dimensional case was studied in 31 ] . The asymptotic behaviour 
of the two-dimensional for some integrals over all complex plane with 
weak singularity was studied in pf . 

A.l Integrals over half-plane with general sta- 
tionary point of fast oscillated exponent 

Here we study an asymptotic behaviour of an integral: 

f f dnAdn . .. 2 

1= —. exp(-z(n 2 + n 2 )), 

J Jq+ I — n 

where |/| — > oo and the domain Vl + = {Re(7 + I — i(l — I)) > 0}. 

Theorem 7 The asymptotic behaviour of the integral I as \p\ — > oo 
has the form: 

I = ^xpH^ + P)) _ 3^ + 0(|r2) l £ n+ . 

I = -™- + 0(|/r 2 ) where l<£9, + . 

The proof. Let us suppose that I G + . Divide the domain £l + 
into three domains: first one is Qf = {Q + \{[\n\ < \l\/2] U [\n — l\ < 
e]}}. This domain hasn't stationary point of phase function of the 
exponent and the singularity of the integrand. Second one is f2^~ = 
{f2 + \[|n| > |/|/2]}. This domain contains the stationary point of the 
phase function. At last, third domain is = {|n — l\ < e}. This 
domain contains the singularity of the integrand. 

Let us integrate by parts over . As a result we obtain: 

dnAdn , , 9 ,. f dn exp(— i(n 2 + n 2 )) 

■exp(-i(n 2 + n 2 )) - 1 



I — n Jan+ I — n —2in 

dn A dn exp(— i(n 2 + n 2 )) 
I — n 2in 2 
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The boundary of the domain includes a large half-circle of radius 
R as R — > oo, a circle of radius e over the point I = n, a half-circle 
of radius |/|/2 and two segments: [i?exp(3i7r/4), \l\ exp(3i7r)/2] and 
[|Z|exp(-wr/4)/2,.Rexp(3z7r)]. 

Consider the integrals over boundary of the domain Qf. The in- 
tegral over large half-circle is equals to zero as R —* oo. The inte- 
gral over half-circle of radius |/|/2 has order (because of oscilla- 
tions and small value of the integrand). The integral over the circle 
at n = I equals to a residue of the integrand multiplied by 2iir as 
e — > 0. The integral over the segments [i?exp(3i7r/4), \l\ exp(3i7r)/2] 
and [|Z| exp(— in /4:)/2, Rexp(3in)] (let us denote their union by L) is 
reduced to the form: 



/ 

J n 



dn exp(3i7r/4) f d\ 



neL -linil - n) —2il J\x\>\i\/2 I - Aexp(i7r/4) 

The second term of (|84j) has to be evaluated as by 0(|Z|~ 3 ). 
Let us consider the integral over Q^. Represent: 

1 I, n s 
= (1 + 



I — n I I — n 
Then the integral over Q^" has to be written as: 

h = / / , = exp(-*(n 2 + n 2 )) = 

J Jttj I — n 

Iff dn/\dn . n _ 2 

■ exp(— i[n +n )) + 



I J Jq+ I — n 

nexp(— i(n 2 + n 2 )). 



Iff dn Adn _ _> > 



I J Jn+ I — n 

Here we integrate by parts the second term. After evaluations we 
obtain: 

I2 = -= / / dn A dn exp(— i(n 2 + n 2 )) + 
I J Jn+ 

1 /•]/|exp(-i7r/4)/2 ^ n 

57/ =^ + o(ir 2 )- 

2Z< 7|/|exp(3i7r/4)/2 t — U 

The integral over O3" as e — > equals to zero. 
Let us sum the obtained asymptotics: 

/ = 2x;r eXp( ~^ 2 + P)) + exp(i7r/2) /"°° dA 

2zZ ' 2iZ 7-oo Zexp(-i7r/4) - A 
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— =. / dnAdnexp(-i(n 2 + n 2 )) + 0(\l\~ 2 ). 
2il Jn+ 

Thus the first statement of the theorem is proved. The second 
statement has to be proved by the same way. 
The theorem is proved. 



A. 2 Asymptotic behaviour of the integral 
with confluent phase function 

In this section we obtain an asymptotic behaviour of an integral as 
\p\ — ► oo 

W+= [ [ 4^£exp(-iw(r)). (85) 
J Jn+ p-r 



Theorem 8 The asymptotic behaviour of the integral (85) where uj(p) 
= 4(p 3 + p 3 ) — v 2 (p + p) as \p\ — > oo and \p ± > has the form: 



£000^ J + p'PoiV" J + Up 1 + 27FZ 15^2 — 



Up 2 



+ 



+0(\p\ +\v\\p\ ), when p€Q + ; 



±^( v 2 ) + ^(v 2 ) 



1 wk+ ^„,2\ TTt 

T2^ 



(86) 



+ 



+0(|p| 3 + |f| 2 |p| 2 ), u>/ien p^^l + . 
Let us prove the theorem. Represent the integral in the form: 



W 



dn A (fnexp(— iu>{n)) + 



+ ^2 



Up 2 



dn A dnexp(— ioj(n))n + 
fi 2 exp(— iui(n)) + 



dn A dn 



+ 



1 



12p 2 



n+ p-n 



n+ p-n 

(Yin 2 — v 2 ) exp(— ioj(n)). 



dn A dn , > 2 



Integrate by parts the last term. As a result we obtain: 
W 



+ 



^o(^) + ^i(^ 2 ) + 



1 



(in (12n 2 — t> 2 ) exp(— iuj(n))) 
12p 2 Jdn+ p — n (— i)(12n 2 — v 2 ) 

o(\ P \- 3 + \v\ 2 \ P r 2 



+ 
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Consider the integral over the boundary of the domain <9f2 + . This 
integral may be considered as a sum of integrals over imaginary axis 
Re(r) = 0, over half-circle of radius R as R — ► oo, Re(r) > and 
over the circle \p — r\ = e as e — > (if p E f2 + ). We can see that the 
integral over the large half-circle tends to zero. As a result we obtain 
the statement of the theorem. 

One more double integral which we need is: 



J Jn+ p-r 

Its asymptotics can be evaluated by the same way as the asymptotics 
of the integral ( |85| ) as \p\ — > oo \p ± v\ > 0: 



jj _ +0{\p\ 2 + \v\ 2 \p\~ 2 ), when p £ 

i<Pto(v 2 )-% + 0(\p\~ 2 + \v\ 2 \p\- 2 ), when p 0+. 

B Reducing of four- multiply integral 
into double integral 



This section is pure technical. Here we show as the four-multiply 
integrals have to be reduced into double integrals over half-plane and 
all complex plane. 

B.l Four-multiply integral with nondegener- 
ate phase of the exponent. 

Let us show that the four-multiply integral may be written as a sum 
of double integrals. Change the variable: n — m = r, then the integral 
J has the form: 





Integrate by part over n. As a result we obtain: 
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dv A dv . , o _o \ \ / r\ • / \ \ 

= exp(«(r +r J) exp(— 2i[rn + rnj) + 

dnn 



— lim f 

Im e^OJu. 



2m e^Oj\i^ n \ =£ Ij - n 

1 f f dr A dr 2 _ 2 

— / / z expmr +r )) exp( — 2i[rn + rn)) — 

m J JC r 



2m 

1 f f dn Adn 



2m J Jc lj — n 



■n x 



— - [ [ dr A dr{— 2i) exp(i(r 2 + r 2 )) exp(— 2i(rn + rn)). 
2m J JC 

First term is equal to zero because of fast oscillating of the integrand, 
second term is equal to residue with sign minus of the integrand at 
n = I, the internal integral in the third term has to be evaluated. 
Finally we obtain: 

exp(i(7 2 + E)) r r dn A dn . .. 2 ^ 

3 = — 2^ / Jc exp( - i( " +fl » - exp(i( '? 

B.2 Four-multiply integral over different half- 
planes with nondegenerate phase function 

Let us consider an integral: 

dn A dfi exp(i(n 2 + n 2 )) 



J_ 



a- 



l — n 2m 

dm A dm 2 -2w 

■ exp(— i{m + m )). 



ln+ n — m 
A partially integrating over n gives: 

I' dn _ exp(z(n 2 + n 2 )) 
+ ^ — n 2i7r 

dm A tim 



n+ n — m 

2 i 



exp(— i(m + m )) 



dn/\dn_ d ex.r>(t(n 2 + n 2 )) 

— ; "-773 — — x 

n- £ — n an 2i7r 

dm A dm . .. 2 -2\<, 
exp(— iym + m J). 



y In 



ln+ n — m 
The derivative is: 

, dn n exp(i(n 2 + n 2 )) 

J u = / — — X 



an- I — n 2m 
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dm A dm 



q+ n — m 

3 f dn 



exp(— i(m 2 + m 2 )) 



Idi 



exp(i(n 2 + n 2 )). 



4 Jan- I - n 

Let us consider integral over part of the boundary in first term : 

,.,Rcx P (3«r/4) dn nexp(i(n 2 + n 2 )) 



X = lim 



Rexp(-i7r/4) I ~ U 

dm A dm 



2itr 



n+ n — m 



exp(— i(m 2 + m 2 )). 



Change the order of integration over m, m and n. The internal integral 
over n has to be evaluate: 



f 

lim / 

R-^oc Jr 



i?exp(3i7r/4) 



dri- 



ll 



cxp(-i7r/4) (I — n)(n — m) 



when I G 



i _|_ m 

-1/2 when / G £1+. 



In last expression of J_ + the integrals over the large half of circle 
|n| = R as i? — ► oo are equal to zero. 

Final formulas for J__i_ have to be written as: 



J- 



-ZTT 



+ il J J n+ dm A dm ^~f™ 2 +™ 2 )) ; W hen Z G 



+ §z7rexp(Y)P + P)) - /exp(i(/ 2 + P))x 



f /„+ dm A dm exp( -^ + " 2)) , when lefl". 



B.3 Four-multiply integral with confluent 
phase function 

Let us reduce an four-multiply integral 



Ji = 



f f dn A dn f f dm A dm 

/ / exp iiuj n) / / - expi— iu(m)) 

J JC p-n J JC n-m 



dm A dm 



p — n J JC n — m 

into sum of double integrals. Denote: 



V(n, v 



2 , exp(iu)(n)) 

2iTT 



//c 



dr A dr 



C n — r 



exp(— iu(r)). 



The function V(p, v 2 ) has to be written as: 
dV 

-— = 12ifi exp(iu>(n))(poo(v 2 ) + 12iexp(iu>(n))(j)oi(v 2 ). 
on 
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We can obtain this formula by changing variable in the integrand 
p = n — r and differentiating the obtained expression of V(n, v 2 ) with 
respect to n and, finally changing the variable back: r = n — p. 
Partial differentiating of the double integral over n and n gives: 

dnAdfi^ r , , f dnnV(n,v 2 ) 

■V(n,v 2 )= hm / — — - + 



C p — n ' R^coJ\ n \ =R p — n 

i dnn Tr/ 2n f f dn A dn _ dV(n, v 2 ) . 
+ hm / V(n,v 2 ) - / / n ^ — -. (87) 



e^o J\ p - n \ =£ p - n ' J JC p-n dn 

Asymptotic behaviour of V(n, v 2 ) as \n\ — ► oo was obtained above. 
Using this asymptotics for the first part of the formula (37) gives 
zero. Second term of the right-hand side of ( |87| ) gives the residue of 
the integrand at p multiplied by 2m. Then we obtain: 

Ji(p,v 2 ) =pV(p,v 2 ) - 4>oo(v 2 ) exp(m>(p)) - 
2 . 1 f f dnAdn. 



-Y2i<j)ax(v )— I I nexp(zu;(n)) + 

2m J JC p-r 

. 2 , 1 f f dn Adn ., . N 
+ iv (poo{v )— / / exp(zcj(n)). 



2m J JC p — r 
On the same way one can evaluate an integral 

I f f dn Adn f f dm A dm 

J, =11 expucjm)) / / - eyp( — m . 

J Jn- p-n J Jn+ n — m 

The difference between the integral J\ and the integral + con- 
sists in the addition terms over boundaries of the domains f2 + and 
in the result. The finally forms are: 
when p 6 

Jf + = -(j)oo(v 2 ) exp(iu)(p)) + 

/* /" dm A dnx 
+pexp(iu;(p)) / / - exp(— iuj(m)) 

J Jn+ p-m 

Iff dn A dn _ 

- / / nexp(iui(n)) + 

2 J 7n- p — n 

2 , , 2 > 1 /" /" dn Adn 



+iv 0J o (u )— / / exp(iw(ra)) 

2«7r 7 7n- p-n 



2 , 1 /" /" dnAdn, 



-12i05i(« )— / / nexp(iw(n)); 

2i7r 7 Vfi- p-n 
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when p fj+: 

T . ,+ / 2\ ft dm A dm . 

J\ = -0oo(v )+P / exp(-w(m)) 
J Jf2+ v + m 



Iff dn Adn _ ... 

n exp(«u;(ra)) 



iw 2 



2 7 in- p - n 
^+ / 2\ 1 /" /" dnAdn ,. , 

^ (w) 2^il-T^T expMn)) 

...... j. 1 /" /" dnf\dn_ 

l2l ^ v y y„_ -^r" expMn)) - 
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